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Abstract 

A general definition of the Spintronics concept of spin-pumping is proposed as generalized forces 
conjugated to the spin degrees of freedom in the framework of the theory of mesoscopic non- 
equilibrium thermodynamics. It is shown that at least three different kinds of spin-pumping 
forces and associated spin-currents can be defined in the most simple spintronics system; the 
Ferromagnetic/Non- Ferromagnetic metal interface. Furthermore, the generalized force associated 
to the ferromagnetic collective variable is also introduced in an equal footing, in order to describe 
the coexistence of the spin of the conduction electrons (paramagnetic spins attached to s-band elec- 
trons) and the ferromagnetic-order parameter. The dynamical coupling between these two kinds 
of magnetic degrees of freedom is presented, and interpreted in terms of spin-transfer effects. 
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I. INTRODUCTION 



Spintronics is a generalization of electronics that takes into account the degrees of free- 
dom related to the spins of the conduction electrons (or the spin of other electric carriers). 
A typical spintronics system is defined by a statistical ensemble of a plurality of electronic 
populations (discriminated by their internal degrees of freedom), that are put out of equi- 
librium in the presence of electric and magnetic forces. The consequence is the creation 
of currents of electric-charge carriers and currents of spins. Spintronics emerged with the 
discovery of giant magnetoresistance in the late 80s l|, and it plays today a crucial role in 
the development of new electronic devices and functionalities. 

The goal of this report is revisiting spintronics on the basis of the theory of Non- 
Equilibrium Thermodynamics {2 6|. The analysis is based on the first and second law 
of thermodynamics, i.e. on the expression of the power dissipated through the different re- 
laxation mechanisms that characterize the system. The description holds at the mesoscopic 
scales, under the hypothesis of local equilibrium extended to internal degrees of freedom 
71, |8|. This work is restricted to classical systems. The extension to quantum systems, and 
especially for the definitions related to "permanent currents" (for which the second law of 
thermodynamics is inoperative) is beyond the scope of this report js]. 



1. Longitudinal spin relaxation 

The typical system to be investigated is a ID wire (with x as coordinate) containing a 
Ferromagnetic/Non- Ferromagnetic metal junction (or a Ferromagnetic/Ferromagnetic metal 
junction). In the most simple cases it can be described by a statistical ensemble of indepen- 
dent electrons that are defined by an effective mass m*, an electric charge e, internal degrees 
of freedom reduced to the spin one half, and a quantization axis (fixed by the magnetization 
of the ferromagnetic layer). The system can then be reduced to two electronic populations 
locally represented by a reservoir of conduction electrons of spin up - defined by a chemical 
potential /it(x) - and a reservoir of conduction electrons of spin down - defined by a chemical 
potential fJ^iix). This is called the two spin-channel model. 

There are then two different ways to put the system out of equilibrium, (i) Applying an 
electric field = ^ ^ will create spin-dependent electric currents J| = cr^E-^ through the 
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wire (where is the conductivity), (ii) applying an effective magnetic field will put the two 
spin populations out-of-equilibrium, creating in turn a flux of spins (noted ip here) in the 
spin configuration space. 

However, the striking point is that the two generalized forces, electric and magnetic, 
are not independent. Applying a voltage difference through the junction generates also a 
stationary flux of spins ip because the electric field difference Ef — E^ = -^^§^ is maintained 
in the junction. The junction produces a non-zero chemical potential difference Afi = fi^—fi^, 
that plays the role of an effective magnetic field, or spin-pumping force conjugated to the flux 
of spins ip. The spin flux can also be written = Sn/r, where r is the spin- flip relaxation 
time and 6n = — is the density of out-of-equilibrium spins. 

As a consequence, applying a voltage difference through the junction generates not only 
stationary spin- dependent electric currents J| but also the stationary spin flux ip. 

If £(5,n-t-,n|) is the energy density of the system (function of the entropy S and the 
density of charge carriers % and n^.), the following canonical definitions holds: T = H, 
/i^ = ^ and A/i = The first relation defines the temperature, the second defines 
the chemical potential, and the last relation defines the pumping force Afx as the chemical 
affinity of the reaction that transforms a conduction electron of spin up into a conduction 
electron of spin down fu, 12|. 

The power dissipated by the system at fixed temperature reads TdS/dt = —J^E^ — Jt^E^ + 
ijj A/j,, where dS/dt is internal entropy production. The relation between the generalized force 
and the generalized flux is imposed by the second law of thermodynamics, dS/dt > 0, and 
is formally expressed by the Onsager relation. 

This model has been first proposed by Johnson and Silsbee [l^, and the introduction of 
;he pumping force Afi in the context of spin-dependent transport is due to Van Son et al. 
1^ . with the description of longitudinal spin-accumulation and giant magnetoresistance at 
a Ferromagnetic/Non- Ferromagnetic interface. This approach was systematized by T. Valet 
and A. Fert in 1993, on the basis of the Boltzmann transport equations fis!]. The model is 
presented below (Section III) in the language of non-equilibrium-thermodynamics. 

Note that in the Spintonics literature, the term spin-current is devoted to the spin- 
dependent electric current 6J'^ = — that flows through the wire, and not to the 
spin- flux ip that flows though the spin-space (e.g. the so-called "Bloch sphere"). 
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2. Spin precession 



However, the description of spin dynamics in the spin-space is not restricted to spin-flip 
relaxation (i.e. "longitudinal" spin relaxation) but it is also characterized by precession 
effects ("transverse spin relaxation") as described e.g. by the Bloch equation of the para- 
magnetic resonance or by the Hanle effect measured in semiconductors [17]. The effect of 
the precession can be taken into account in terms of diffusion of the transverse components 



of the spin. This generalization was called transverse spin-accumulation [18|, ll9|], and was 
introduced much later in the context of spin-transfer experiments (see below). 

As a consequence, there is another way to drive the spins out-of-equilibrium with the 
use of the transverse spin pumping force Afij_ that generates the transverse spin flux ijj^. 
Also in that case, the application of a potential difference trough the junction gives rise to 
non-zero A/i_|_, and the transverse spin-polarized current 6J± = cr±^^^ is produced at the 
interface. The tranverse spin flux can also be expressed in terms of transverse relaxation 
time r_L with 'iIj± = 6n±/T±. These recent developments are presented in Section IV on equal 
footing with the longitudinal spin - flip. 

3. Band structure and s-d relaxation 

Furthermore, beyond the spin internal degrees of freedom, it is important to push the 
description to a more realistic situation that takes into account the ferromagnetic specificity 
of the material (and not only the paramagnetic properties of the conduction electrons). 
Indeed, the out-of-equilibrium magnetization described within the two-channel model is 
given by 6m = gfisi^^ ~ "^4.) (where fiB is the Bohr magneton and g the Lande factor). This 
system is paramagnetic. The contribution Sm to the total magnetization is superimposed 
';o the ferromagnetic collective variable M (which is essentially due to the d band electrons 
20l|). In terms of transport properties, the quasi particles s and d are defined by effective 
masses m* and m^, i.e. by supplementary internal degrees of freedom that takes into account 
the coupling to the periodic lattice. 

In line with the pioneering works of Mott [l6^, we will consider a simple generalization 
of the two channel model that takes into account the ferromagnetic nature of the 3d metals 
with enlarging the internal degrees of freedom to four electronic populations: the conduction 
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electrons of the s band for up and down spins, and the conduction electrons of the d band 
for up and down spins. We have then a four-channel model, in which two kinds of spin flux 
ip and ipsd should be defined for both interband and intraband spin- dependent relaxation. A 
third kind of pumping force can then be introduced with the interband chemical potential 
difference A/x^^ j21|. This is performed in Section V. 



4- Ferromagnetic collective variable 

Nevertheless, there is something missing in the above description of a ferromagnetic 
junction: the ferromagnetic collective variable M has not been introduced explicitly. Ac- 
cordingly, the next section below (Section II) presents a derivation of the equation of the dy- 
namics for the ferromagnetic variable (i.e. the Landau-Lifshitz-Gilbert equation) performed 
with the introduction of a generalized force Ve/i^ = —H^ff thermodynamically conjugated 
to the ferromagnetic degrees of freedom (/i^ is the ferromagnetic chemical potential and 
S is the magnetic configuration space). The effective magnetic field H^ff can also legit- 
imately pretend to the appellation "spin-pumping force". This generalized ferromagnetic 
force generates the current in the magnetic configuration space. 

Yet, many experiments have shown that it is possible to switch the magnetization 
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25| or to generate ferromagnetic entropy 26| while injecting spin-polarized currents in 



Ferromagnetic/Non- Ferromagnetic junctions. The corresponding effect - called spin-transfer 
- shows that the spin-dependent electronic transport coefficients are coupled to the transport 
coefficient of the ferromagnet. In other terms, the ferromagnetic current generated by H^ff 
and the spin-polarized current generated by A/i, Afi±, or Aftgd are coupled. 

The dynamical coupling that occurs between the current of spins and the current of 
ferromagnetic moments is discussed in Sec. VI. It leads to define a spin-transfer effect for all 
;he spin-pumping sources we have previousely identified : longitudinal Q, Q , transverse 
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29| . and s-d interband relaxation 30|, |31| 



II. INTRODUCTION OF THE FERROMAGNETIC DEGREES OF FREEDOM 

In this section, we focus on a uniform ferromagnetic moment M = MgUr defined with 
radial unit vector Ur and the magnetization at saturation Mg. It is called macrospin, in 
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FIG. 1: Illustration of the S configuration space of the ferromagnetic variable. The two potential 
minima and the ferromagnetic currents are sketched 



opposition to the microscopic spins attached to atonies or electrons. In order to treat 
statistically the time dependence of ferromagnetic degrees of freedom m contacted to a heat 
bath, the ergodic property is used. It allows working with a statistical ensemble of a large 
number of ferromagnetic moments m±dm that defines a surface {6^6, (p±dip} on the sphere 
S of radius M^. The corresponding de nsity p^{6,ip) is then identified with the statistical 



distribution of ferromagnetic moments [32|. The mean value of the magnetization is M. 
The introduction of the density is justified by the nanoscopic size of the magnetic single 
domain, for which the fluctuations play a major role. To that point of view, the system is 



mesoscopic. Accordingly 33|, the ferromagnetic chemical potential takes the general form 
/i^ = kTln{p^) + where is the usual ferromagnetic potential (deduced e.g. from the 
quasi static hysteresis loops) and the first term accounts for diffusion. On the other hand 
the current of ferromagnetic moments, = p^dur/dt, is confined on the surface of the 
sphere. 

The power dissipated by the ferromagnetic system is given by the corresponding inter- 
nal entropy production and is formed by the product of the generalized flux by the 
generalized force. Assuming a uniform temperature T we have: 
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The application of the second law of thermodynamics dS^ /dt > allows the transport 
equation to be deduced by writing the relation that links the generalized flux (the current 
J^) of the extensive variables under consideration and to the generalized force defined in the 
corresponding space S. Both quantities, flux and forces, are related by the Onsager matrix 
of the transport coefficients C: 



.J^ = -CV^li^ (2) 

This is the simplest form of the well-known Landau- Lifshitz equation (see below). We 
started from the hypothesis that the magnetic domain is uniform: the modulus of the 
magnetization is conserved. The trajectory of the magnetization is then confined on the 
surface of a sphere of radius Ms, and the fiow is a two component vector defined with the 
unit vectors {u^jUe} of E. Accordingly, the Onsager matrix is a 2x2 matrix defined by four 
transport coefficients {L00, Lo^p, L^po, L^^p] . Furthermore, the Onsager reciprocity relations 
impose that Lg^p = —L^g. Assuming that the dissipation is isotropic, we have Log = L<^^. 
Let us now introduce a dimensionless supplementary coefficient a, which is the ratio of the 
off-diagonal to the diagonal coefficients: a = Lg^p/Lgg. In conclusion, the ferromagnetic 
kinetic equation is defined by two ferromagnetic transport coefficients Lg^p = p^Lp and a: 




£ = p'Lp\ ^ I (3) 

On the other hand, the generalized force Vs/U^, thermodynamically conjugated to the 
magnetization, is the effective magnetic field H^ff = — VsM^- It is a generalization in the 
sense that this effective field includes the diffusive term [SJ] that was first introduced by 
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Brown in the rotational Fokker-Planck equation 

Actually, it could be rather surprising to claim that Eq. (jj]) is the "well-known LL 
equation" . However, it is sufficient to rewrite Eq. ([2]) in 3D space with re-introducing the 
radial unit vector Ur = (1,0,0) of the reference frame {ur,ug,u,p}, and recalling that the 
current is the density multiplied by the velocity = p^dur/dt, to recover the traditional 
LL equation from Eq. ([2]) and Eq. ([3]): 



dUf. 
~dt 



-Lp jwr X Heff + aUj. X X -ffe//| j (4) 
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Furthermore, it is well-known that LL equation is equivalent to the phenomenological 



Gilbert 
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equation, that defines the magnetic damping coefficient r]: 

— ^,M.\^H,„-„—'j (5) 

where 7 is the gyromagnetic ratio. The equivalence between the two equations defines the 
coefficients a and L^? as a function of the coefficients 1] and 7: a = rj'^Ms is the dimensionless 
damping coefficent and Li? is defined by the relation: 

^ (6) 



" Ms (1 + a^) 

The above approach can be applied to microscopic spins (e.g. for the derivation of the 
Bloch equation), but it should be generalized to the case in which the modulus of the 
magnetization is not constant (3X3 matrix) and the damping is not necessarily isotrope 
{L0 7^ L^^^). This task is beyond the scope of this report. 

III. TWO SPIN-CHANNEL MODEL 

In this section, we only focus on the spin-dependent electric transport only, and on the 
two-channel model sketched in the introduction. The corresponding electric wire is defined 
along the x axis, with a section unity: the relevant configuration space is the one- dimensional 
real space M. 

The conservation laws write: 



IT ~ ~~dx ^ /n\ 

dt ~ dx ^ ^ 



where r2-|- and are the densities of charge carriers in the channels {t;i}; and the spin- 
dependent relaxation is taken into account by the fiux ip. This is the velocity of the reaction 
(or relaxation of the spin-dependent internal variable) that transforms a conduction electron 
'I" into a conduction electron J,. This relaxation is formally equivalent to a chemical reaction. 



driven by the chemical affinity A/x = /i^ — [HI]. The power dissipated by the system then 
reads: 
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T 



^ e dx 



(8) 



dt ^ e dx 

The corresponding kinetic equations are deduced from the second law of thermodynamics 
after introducing a supplementary Onsager coefficient L. 



X 



cr-[- dfl^ 
e dx 



(9) 



^ e dx 
tp = LAfi 

The set of equations Eqs @ is sufficient and necessary in order to describe, in the 
stationary regime, spin-accumulation effects and any non-equilibrium contribution to the 
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2l|. 



resistance due to relaxation (A/i 7^ 0) occurring at an interface 

Eq. on]) shows that the spin-dependent electric currents J| and the spin flux ip are not 
independent. Accordingly, it is more convenient to rewrite Eq. ([9]) as a function of the 
variables Afi and fiQ = fi^ + fi^. Let us deflne the conductivity asymmetry by the parameter 
/3 such that /3 = -^^^ and the mean conductivity 2cro = cr^ + cr^. On the other hand, 
the spin-polarized electric current is 6J^ = J| — J| and the spin- independent current is 
Jq = J| + Jf . In this new system of equations, the Onsager matrix re- writes: 



^0 



( 



\ ^ ) 



\ ( ~1 d^P-" \ 
e dx 

-1 dji-l 
e dx 



\ 



J 



(10) 



ao /3o-o 
/3(Jo do 
L 

The current 5J^ is called "spin current" or "pure spin-current" in the spintronics litera- 
ture. 

The system of equations Eq. (ITOi) allows the diffusion equation for A/i(x) to be derived 



for the stationary conditions ^ = and 



dx 



dx 



dx'^ 



''diff 



(11) 



where ij^^ 



2eL 



<To(l-/9")- 



And the non-equilibrium magnetoresistance produced by the 



interface writes (see the details of the derivation in references 
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FIG. 2: Spin-accumulation An{x) trough a Non-Ferromagnetic/Ferromagnetic/Non-Ferromagntic 
junction with typical size of the order of the spin diffusion length 

where the measurement points A and B are located far enough in each side of the interface 
so that Ajj^A) = A/i^(5) = 0. 

The spin diffusion length is typically some tens of nanometers in ferromagnetic metals, so 
that the non-equilibrium magnetoresistance requires thin films (current-in-plane geometry) 
or nanostructured pillars (current-perpendicular-to-the-plane geometry) to be exploited. 



IV. SPIN PRECESSION 



The description proposed above with a spin-dependent internal variable that takes the 
two spin values {t, 1} is not able to take into account the precession of the spins occurring in 
a magnetic field, and observed with electronic resonance or Hanle effects. In the case of the 
processes that lead to spin accumulation and giant magnetoresistance, the mean values are 
averaged out over the spin-diffusion length, so that the precession of the spin is not relevant. 
However, this is no longer the case in a quasi-ballistic regime close enough to the interface. 

In order to take into account these quasi-ballistic effects (i.e. sub-nanometric scales in 
metalic devices), the two-channel model has been recently generalized to transverse spin- 



accumulation in the context of spin-transfer-torque investigations 
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37( 1 . The transverse 



spin-accumulation is introduced with the corresponding current 6J1 and the corresponding 
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chemical potential A/x^. Transverse means here that the spin density is considered in the 
plan perpendicular to the quantification axis | that defines the spin up and spin down in 
the two-channel model. 

The conservation laws writes: 



d6n± 



(13) 



dt dx 

The transverse spin fiux iIj± can be expressed with a transverse relaxation time t±: ip± = 

where n_|_ is the density of transverse spins. 
Note that the two potentials A/i*^ and A/x^ are defined at very different length scales and it 
is necessary to refer to quantur a app roaches in order to understand the physical signification 



of the transverse parameters 29 
dissipated power is 



38| . The corresponding transverse contribution to the 



T 



-6J' 



^ + ^xA/ij 



dt edx 

Putting all together, we have the following Onsager relations for the electric system: 



(14) 
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(15) 



V. THE ROLE OF THE d | ELECTRONIC SUBBAND 



A justification of the spin-dependent conductivity asymmetry /3 7^ in ferromagnetic 



metals has been proposed by N. Mott in 1936 [16|], on the bases of the newly discovered 
band-structure approach. In the Mott description, the observed transport properties (e.g. 
the huge resistivity of Ni below the Curie temperature) have been accounted for by the 
existence of four electronic populations: the conduction electrons of spin up and down (s ^ 
and s D of the s band and the conduction electrons of spin up and down {d ^ and d l) of 
the d band. The argument is based on the fact that the contribution to the resistivity of 
the s — d interband scattering is higher than the contribution of the s intraband scattering. 
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In the ferromagnetic 3d metal, the d t band is full so that the relaxation channel of s t 
electrons to f band is blocked (according to the Fermi golden rule, the relaxation rate 
is proportional to the density of states in the final d J, band). Furthermore, the spin- flip 
interband relaxation is too energetic to be efficient (the relaxation s f to 4 is negligible). 
As a consequence, the s l are more scattered that the s t, and the conductivities of the two 
channels is asymmetric: cr^ > cr^.. This mechanism is also responsible for the anisotropic 
magnetoresistance [39^. The necessity of enlarging the internal degrees of freedom to the 
band structure leads to enrich the concept of spin-currents and the spin-pumping force. 

Using the notations introduced in the previous sections, the total current Jt is composed 
by the three currents for each channel : Jt = Js-\ + ^^4. + Jdi {Jd-t = because the band is 
full). The relaxation rate tpsd is introduced to account for s — d spin-conserved scattering, 
and the relaxation rate ips, is introduced in order to account for previously defined spin-flip 
scattering. Assuming that all channels are in steady states, the conservation law write: 



drit 


dJt 


dt 


dx 




dJs^ 


dt 


dx 




dJsi 


dt 


dx 


dudi 


dJdi ^ 


at 


dx 







i^sd + i^s 
ipsd = 



(16) 



where nt,ns^,ns^,ndi are respectively the total densities of particles and the density of 
particles in the s fj 4, J, channels. The conjugate (intensive) variables are the chemical 
potentials {fig^, fJ'si, /^df f^di}- The application of the first and second laws of thermodynam- 
ics allows us to deduce the Onsager relations of the system : 



J si 
Jdi 

i^sd 



CTsi dfi. 



si 



e dx 
e dx 

CTdi dfidi 



(17) 



e dx 

Lsd ifJ'si — fJ'di) 
Ls ifJ'st — fJ'si) 

where the conductivity of each channel {cTs^, <Jsi, o'dt, o'di} has been introduced. The 
first four equations are the Ohm's law applied to each channel, and the two last equations 
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introduce new Onsager transport coefficients, Lsdi and Lg, that respectively describe the 
s — d relaxation (I) for minority spins under the action of the chemical potential difference 
A/i| = /is4./2 — fidi and the spin-flip relaxation (II) under spin pumping A/ig = fi.^ — Hsi/'^- 
The Onsager coefficients are proportional to the corresponding relaxation times [12|. 

In the same manner as performed in Section III, the equations of conservation Eqs. f ll6p 
and the Onsager equations Eqs. ( IT7I) lead to the two coupled diffusion equations : 

f S^A/i, 1 1 
d'^lis 1 , 1 , ^ ^ 

^•^ ^sd ''sf 

where the four diffusion lengths {Igd, As, hf, ^sd} are given as a function of the transport 
coefficients in reference jl2 |. 

This three-channel model brings to light the interplay between band mismatch effects and 
spin accumulation, in a diffusive approach. The resolution of the coupled diffusion equations 



is discussed elsewhere |l2 |. 



VI. DERIVATION OF SPIN TRANSFER DUE TO SPIN-PUMPING FORCES 



In usual experimental configurations for spin-transfer, an electric current is injected in 



a ferromagnet through an interface (in series or in non-local configuration [40|]) and the 
magnetoresistance, i.e. the potential drop Eq. (fT2l) allows the magnetization states to be 
measured. The effect of strong electric currents on the magnetization states can then be 
observed. In such a configuration, the two sub-systems S and the spin-polarized current 
described in the previous sections exchange magnetic moments at the junctions and both 
are open systems. 

In order to describe the dynamics of the ferromagnetic degrees of freedom, we have to 
deal with a closed system. The system of interest is now the ferromagnetic system that 
includes spin-accumulation effects at the junctions. 

For the sake of simphcity, we treat in the following a unique decoupled spin-dependent 
process A/i and ip that includes s — d relaxation. 

This total ferromagnetic system is such that the density of ferromagnetic moments and 
the total ferromagnetic flux J^^ are related by the conservation law: dpf^^/dt = —div^J^^. 
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The initial configuration space of magnetic moments is tlien extended to ID real space 
parametrized by the internal variable S The important point here is that the internal 

variable is spin dependent, and related to the ferromagnetic space S (e.g. through spin-flip 
or s — d relaxation and the corresponding spin accumulation). This accounts for the coupling, 
i.e. the transfer, of magnetic moments between the two sub-systems. 

The dissipation is given by the internal power dissipated in the total system T dS / dt : 



+ + (19) 



edx edx edx 

Where the first term in the right hand side is the power dissipated by the total ferro- 
magnetic sub-system (including the ferromagnetic contribution due to spin-transfer), the 
two following terms are the power dissipated by spin-dependent electric transport, and the 
fourth term is the spin-independent Joule heating. The last term is the power dissipated by 
spin-flip or s — d relaxation. 

In Eq. ( fT9l) . the vectors are defined on the sphere S with the help of two angles 6 and 



(p. The total ferromagnetic current jf^^ = itoti.'^v jtot//^s includes the contribution due 
to spin-accumulation mechanisms. The chemical potential /x^ accounts for the energy of a 
ferromagnetic layer. On the other hand, the system is contacted to electric reservoirs with 
the electric currents and the corresponding chemical potentials. Applying the second law of 
thermodynamics, we obtain the general Onsager relations: 



' J toil. > 
?^ 

J tot II 

\ ^ J 



—pLp pLpa /// 

l± ax 

/// ao 

Lx 

L/ 



-1 dAf^l 
e dx 
-1 dAii-' 
e dx 



(20) 



V A/i / 

All coefficients were defined in the previous sections, except the two cross-coefficients 
{/_!_,///}, introduced in this model as spin-transfer coefficients. The coefficients are 
given by the Onsager reciprocity relations. We assumed that the other cross-coefficents are 
zero or negligible. 

The total ferromagnetic current can be written after integrating over the volume v of the 
ferromagnetic layer of section unity and the spin accumulation zone. This volume is such 
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that v = dz, where x = A and x = B are two sections close to the interface but far 
enough with respect to the diffusion lengths. We assume here that the diffusion lengths are 
much smaller than the width of the ferromagnetic layer in order to simplify the calculation: 
the volume of the ferromagnet is identified as v. Let us define X as the correction due to 
the spin-transfer after integrating over the volume v deduced from the two first equations 
of the matrix equation Eq. (l20l) .: 

jl, = CvV,:^i^ + pvX (21) 
where C is the matrix defined in Eq. 

The assumption of constant modulus of the magnetization imposes that X is confined on 
the surface of the sphere S. The Helmoltz decomposition theorem can then be applied: the 
vector X can be decomposed in a unique way with the introduction of two potentials x ^^^^ 
$ (i.e. a potential vector) such that: 

X = X Vs$ + VsX (22) 

where the first term is divergenceless and the second term is curlless (i. e. non conserva- 
tive). The total correction to the Landau-Lifshitz-Gilbert equation writes 

Jit = vCVj^fi^ + Vex + Mr X Ve$ (23) 
The generalized LLG takes the form: 

where Heff = — V/i^ is the usual effective ferromagnetic field. 

Eq. is a generalized LLG equation that includes the effect of spin-pumping Afi and 
Afi± (the contribution of Afisd has not been added to the Onsager relations Eq. (!20!) for the 
sake of simplicity). Qualitiatively, the most important point is that the dynamics should 
now be described by the introduction of two potentials, or two magnetic fields H^ff — ^f^ for 
the precession and i?e// — ^0 for the longitudinal relaxation, instead of a single field H^ff 
in the usual case. Both fields appearing in the LLG equation contains a current dependence. 

The corresponding Fokker-Planck equation is obtained by inserting the expression of Jtot 
into the conservation equation: dpf^^/dt = —div-^J^^. 
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VII. CONCLUSION 



The basic concepts of Spintronics - spin-pumping force and spin currents - have been 
presented on the basis of the theory of Non-Equihbrium Thermodynamics. It has been 
shown that different relaxation mechanisms can be invoked, each of which defines a specific 
spin-pumping force and a specific spin flux. Four mechanisms have been investigated within 
this formalism: spin-flip scattering, spin precession, s — d scattering, and ferromagnetic 
relaxation (with both longitudinal relaxation and precessional motion). 

This approach shows that the apphcation of a voltage difference through a 
Ferromagnetic/Non- Ferromagnetic junction leads to the creation of spin-pumping and spin- 
flux, which in turn leads to non-equilibrium interface resistance or spin-transfer to the fer- 
romagnetic collective variable. The derivation proposed for the spin- accumulation formulae, 
dynamics of the magnetization, and spin-transfer effects, are based on the expression of the 
entropy production and the second law of thermodynamics. In principle, it is possible to 
obtain the same results in a more elegant way, by using the Prigogines theorem on minimal 
entropy production. This work rest to be performed. 
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